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ABSTRACT 


A  theoretical  and  experimental  investigation  is 
presented  for  unsteady  free  convection  adjacent  to  a 
heated  horizontal  plate  subject  to  temporal  variations 
of  the  surface  temperature  gradient.  The  equations  of 
motion  describing  the  flow  near  the  surface  were  demon¬ 
strated  to  be  of  a  linear  nature  provided  inertia  and 
advective  effects  could  be  neglected.  On  the  basis 
of  tests  conducted  in  water,  glycerine  and  air  these 
approximations  appear  to  be  valid  over  a  wide  range  of 


Prandtl  numbers . 
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NOMENCLATURE 


x,X 

Y,Y 

u,U 


v ,  V 


F,it 

T  ,t 

P/P/Q 

o 

p 

pd 

4>,*  ,e 

V 


/C 

P 

P 

c 

p 

g 

M- 

a 

Os 

G 

H 

c 


distance  parallel  to  surface 
distance  normal  to  surface 
velocity  in  X  direction 
velocity  in  Y  direction 
stream  function 
time 

pressure 

hydrostatic  pressure 

pressure  due  to  viscous  or  thermal  effects 

temperature 

momentum  diffusivity 

thermal  diffusivity 

thermal  expansion  coefficient 

density 

specific  heat  at  constant  pressure 

gravitational  acceleration 

absolute  viscosity 

Prandtl  number 

Ostrach  number 

body  force  per  unit  mass 

dissipation  function 

constant  with  dimension  of  (temperature)/ (distance) 
(time) 
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In  i  v'  .  .  .  ' 


ix 


n 


€ '  C 


a, A 


^2 
CD, 7 

b 

z 

r 


.  m  _ 

1  erfc 


s 

A 

V 


0 


d 

o 

w 

CO 


surface  temperature  index 
unknown  functions 
arbitrary  constants 
-«4/2  ccx  =  n 

roots  of  indicial  equations 
coefficient  in  Maclaurin  expansion 
arbitrary  function 
gamma  function 

m'th  integral  of  the  complementary  error  function 

r 

Laplace  transform  variable 
finite  change 
Laplacian  operator 


denotes  differentiation  with  respect  to  rj 


SUBSCRIPTS 


departure 

characteristic 

surface 

infinity 


Other  symbols  are  defined  where  they  appear*, 
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CHAPTER  I 


INTRODUCTION 


If  a  shallow  horizontal  layer  of  fluid  is  heated  iso- 

thermally  from  below,  there  is  no  immediate  visible  motion 

of  the  fluid  and,  initially,  heat  is  transferred  from  the 

lower  region  to  the  upper  region  by  conduction  alone.  This 

can  continue  if  the  temperature  differences  in  the  fluid 

are  small,  but  once  a  certain  vertical  temperature  gradient 

★ 

is  exceeded  cellular  motion  occurs.  This  problem  of  the 
motion  of  a  fluid  heated  from  below  in  a  uniform  body  force 
field  has  been  the  subject  of  considerable  study  since  the 
pioneering  work  of  Lord  Rayleigh  (1)  in  relation  to  ex¬ 
periments  made  by  Benard  (2)  in  1900.  He  and  many  authors 
since  then  (3)  have  considered  the  problem  as  one  of  dynamic 
stability.  The  onset  of  instability  is  determined  by  a 
parameter,  the  Rayleigh  number,  which  is  independent  of 
the  cell  shape. 

Flow  may  also  be  established  above  and  below  heated 
isothermal  surfaces  by  a  thermal  buoyancy  effect  (promoted 
by  the  existence  of  body  forces)  due  to  the  presence  of 
leading  edges.  These  edges  set  up  horizontal  temperature 


*  It  is  interesting  to  note  that  a  like  cellular  structure 
may  be  obtained  in  soapy  water  cooling  from  the  surface. 

**  Numbers  in  brackets  denote  references  given  on  page  67  . 
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gradients  within  the  fluid  which  cause  corresponding  density 
gradients  and  hence  pressure  variations  such  that  fluid  flow 
results  and  a  transfer  of  heat  occurs.  Stewartson  (4)  has 
analyzed  the  problem  of  steady,  laminar,  natural  flow  sur¬ 
rounding  such  isothermal  plates  by  a  similarity  technique. 

A  solution  was  obtained  for  heated  plates  facing  downwards 
or  cooled  plates  facing  upwards,  but  in  contrast,  for  heated 
plates  facing  upwards  and  cooled  plates  facing  downwards  a 
mathematical  inconsistency  was  obtained  due  to  an  improperly 
posed  problem.  Levy  (5)  has  tackled  a  similar  problem  using 
the  Von  Karman  momentum  integral  method.  Sugawara  and 
Michiyochi  (6)  used  an  identical  approach  to  calculate  the 
distribution  of  the  local  heat  transfer  coefficient  on  both 
the  upper  and  lower  surfaces  of  a  very  thin  ellipse.  Their 
distribution  agrees  very  well  with  the  experimental  work  of 
Weise  (7)  who  determined  local  heat  transfer  coefficients  on 
rectangular  plates  in  air.  Weise  was  also  the  first  to  show 
that  flow  does  exist  on  horizontal  plates ,  by  means  of 
schlieren  photography. 

On  a  larger  scale,  in  the  field  of  micrometeorology, 
extensive  theoretical  and  experimental  work  on  horizontal 
surfaces  has  also  been  undertaken  with  specific  reference  to 
unsteady  state  conditions  and  turbulent  flow.  Some  of  the 
problems  of  interest  have  been  the  generation  of  mountain 
and  valley  winds  (i.e.  antitriptic  winds)  as  analyzed  by 


V 

_ 
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Jeffreys  (8)  and  the  underlying  principles  of  the  sea  breeze 
as  considered  by  Pearce,  Estoque  and  Fisher  (9,10,11,12). 

In  the  absence  of  a  leading  edge  an  isothermal  surface 
could  not  generate  flow.  But  if  such  a  plate  is  heated  dif¬ 
ferentially,  the  imposed  temperature  gradients  again  produce 
density  gradients  and  hence  pressure  variations,  resulting 
in  motion  adjacent  to  the  surface.  Here  it  is  noticed  that 
a  uniform  body  force  field  is  acting  perpendicular  to  the 
surface  but  there  is  no  particular  restriction  on  heating  or 
cooling  or  on  whether  the  surface  faces  upwards  or  downwards. 
Assuming  thermal  (and  hydrodynamic)  stability,  it  is  evident 
that  laminar  flow  will  exist  if  temperature  gradients  are 
moderate  for  either  steady  or  nonsteady  state  conditions. 

In  the  investigation  to  follow,  the  problem  of  laminar 
natural  convection  adjacent  to  a  large  horizontal  surface 
is  considered.  In  particular,  the  flow  adjacent  to  a  plate 
whose  surface  temperature  gradient  is  a  simple  power  law 
function  of  time  is  examined  both  theoretically  and  experi¬ 
mentally.  The  purpose  of  the  thesis  is:  first,  to  determine 
temperature,  pressure  and  velocity  expressions  describing 
the  flow  on  surfaces  without  leading  edges?  second,  to  veri¬ 
fy  the  validity  of  approximations  made  in  the  theoretical 
analysis  (by  performing  experiments  in  water,  glycerine  and 
air) . 


.  ....  i  *3„x.  ;  n  a  }i  3»S  - ■  »f  ;  <*l  '  a  b^°° 

* 

.  *  ■  ;•  .,  -  ■  =  i  ™  6 


CHAPTER  II 


THEORETICAL  ANALYSIS 


2 . 1  GOVERNING  EQUATIONS 


The  problem  of  free  convection  from  an  infinitely 
large  horizontal  plate  is  considered.  An  attempt  is  made 
to  reduce  the  Navier-Stokes  equation  and  the  corresponding 
energy  equation  to  linear  forms. 

2.1-1  Formulation  of  Problem 


Laminar  incompressible  flows  due  to  heat  conduction  from 
a  surface  to  a  fluid  in  a  uniform  body  force  field  are  des¬ 
cribed  by  the  following  equations? 


Mass  ? 


0 


* 


Momentum  ? 

DV 

Dt 


+  G 


+ 


2~ 
W  V 


2.1-1 


Energy? 


D0 

Dt 


KV2©  +  H 

P 


*  Ignoring  an  extremely  small  contribution  due  to  the 


material  derivative. 


' 


i  r  f. 
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Consider  the  application  of  these  equations  to  the 
two-dimensional  semi-infinite  region  illustrated  in  figure 
2.1.  If  temperature  changes  are  moderate,  all  properties 
except  density  may  be  taken  as  constant  (13) »  The  import¬ 
ance  of  a  variable  density  is  that  it  permits  motion  promoted 


FIG.  2.1  COORDINATE  SYSTEM 


by  buoyancy  forces  in  the  gravitational  field  of  the  earth. 
Density  is  taken  as  a  function  of  temperature  only,  thereby 
permitting  (3,  the  coefficient  of  expansion,  to  be  defined 
from  the  equation  of  state 


p  -  p (9)  , 


in  accordance  with  the  application  of  the  Boussinesqian 
approximations 


oo 


and 


p  -  p 


oo 


^  0 


Thermal  radiation  may  be  present  but  because  most  fluids  have 
a  low  transmittance  (except  gases) ,  its  contribution  is  taken 
to  be  negligible. 


■ 
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Considering  the  surface  temperature  as  a  function  of  X 
and  time,  one  -would  expect  a  boundary  layer  to  form.  The 
resulting  change  in  pressure  will  be  a  departure  from  the 
hydrostatic  pressure  depending  on  whether  the  plate  is  heated 
or  cooled.  Taking 

P  =  P°  +  P, 

d 

permits  equations  2.1-1  to  be  reduced  to  the  following  two- 
dimensional  forms; 


Assuming  no  slip  at  the  surface,  and  since  flow  does 
not  occur  outside  the  boundary  layer,  the  problem  is  sub¬ 
ject  to  the  following  initial  and  boundary  conditions s 


, 


7 


t  *  Os 

(Y 

Y  -  0; 

(t 

Y  —  co  s 

(t 

>  o) 

u  = 

>0) 

u  = 

>0) 

u  = 

9  ~  0 

0,  0  -  0w  (X,t) 

9  "  ?d  ~  0 

2.1-3 


2.1-2  Normalization  of  Equations 


Traditional  methods  failed  to  normalize  equations 

2.1-2  for  these  required  more  than  one  known  characteristic 
quantity.  The  only  characteristic  parameter  available  in 
this  problem  is  the  rate  of  change  of  surface  tempera¬ 
ture  gradient  with  time  (c) .  A  procedure  used  by  Heliums 
and  Churchill  (14)  which  defines  characteristic  quantities 
in  terms  of  fluid  properties  does  permit  the  normalization 
of  equations  2.1-2. 

If  the  dimensional  quantities  are  written  as 

U  =  uU  ,  V  =  W  ,  t  -  Tt  ,  X  =  xX 
o  o  o  o 

2.1-4 

Y  -  yY  ,  P,  —  pP ,  and  9  =  4>ct  X  , 

2  o  d  c  do  o  o 


and  substituted  into  equations  2.1-2,  these  equations  then 


become 


<3u  /XoVo\  <3v 

SS  +  Iyu-J  S7 

x  O  oy  2 


0 


* 


Taking  n  =  1  for  simplicity. 


t  b  ■  ■  1 

* 

•» 


c)u 

5^ 


8 


/toPdoN\ 

vpvoV 


2.1-5 


+ 


gct2X 
o  o 


<D 


d2<t> 

dx2 


+ 


+ 


2 


C  cX  Y 
p  o 


( 


vu 


c  cX  Y 
p  o  o 


2 


VU  V 
o  o 

2 

cX  x 
P  o 


P  o 


Assuming,  for  completeness,  that  all  coefficients  have 

2 

the  same  order  of  magnitude,  each  group  except  /ctQ/Xo  and 


n_&' 


... 


9 


2 

/ctQ/Yo  may  be  equated  to  unity.  This  will  result  in  a  system 
of  algebraic  equations,  the  solutions  of  which  permit  ref¬ 
erence  quantities  to  be  defined  in  terms  of  fluid  properties , 

2  2 

and  the  characteristic  quantity  (c)  .  /ct q/Xq  and  *ct^Y 

2  2 

are  excluded  because  they  combine  with  vt  /X  and  vt  /Y 

CT  O  O  O 

to  give  «/v  =  1  which  contains  no  reference  quantity. 

By  not  forcing  these  terms  to  unity,  they  appear  as  a  dim¬ 
ensionless  group  (a)  in  the  boundary  layer  equations.  Solving 
the  above  mentioned  equations  results  in 


x0  =  Yo  =  <vto)1/2  '  Uo  -  Vo  =  (v/to)1/2 

pdo  “  M/t0  /  a"d  to  =  Ogc)  “  1/3 


2.1-6 


Utilizing  these  definitions ,  it  may  be  shown  that  the  last 
five  terms  of  the  energy  equations  again  combine  into  an 
expression  which  contains  no  reference  quantity.  Hence  sub¬ 
stituting  2.1-6  into  2.1-4  gives 


x  * 


X, 


u 


(Pgcv3)1//6  ' 


V 


Ogcv3)V’6 


(pgo)1/3  t. 


(|x3Pgo) 1/3 


0  . 


and  $ 


'  :  i 

*;•  r.:  ■  L  '  *  '  •:  itj  f  :a  •<*  ©*  a-  AS;  /»  *P  :-v  ;v.:- 

> 
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As  a  consequence,  the  problem  may  be  rephrased  to  one  of 
obtaining  a  solution  to  the  following  normalized  equations s 


c)U  ,  bv 
dx  by 


0 


c)u 


+ 


bv 


+ 


<t> 


a 

by 


b± 

c3t 


+ 


1  fb\_  b^± 

°  W  by2 


2,1-7 


+ 


+ 


fbu  bv\ 
\by  bxj 


where  Os 


vV>g)5/6 

71767 


Initial  and  boundary  conditions  now  become  s 


*  This  ratio  is  analogous  to  the  Eckert  number  in  forced 

convection.  It  is  suggested  that  it  be  called  the  Ostrach 


Number,  Os. 


' 


<<  r*.i ic  .w  cbn 
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T  =  Os 

(y>o) 

u  “  <J)  -  0 

y  =  Os 

(t>  0) 

u  ss  o,  <t>  -  <t>  (x,t)  2.1-8 

w 

y  =  oo : 

(t>  0) 

U  ss  (J)  SS  p  ss  0 

2.1-3  Order  of  Magnitude  Analysis 


It  has  been  observed  for  boundary  layer  flows  that 

Y  «X. 

Since  X  =  Y  ,  it  follows  that 
o  o 

y<<x, 

and  hence  from  the  continuity  equation  2.1-7 

v  «u. 


Although  an  order  of  magnitude  analysis  may  be  made  on 

* 

the  individual  equations  2.1-7  ,  solutions  of  the  combined 
momentum  equation  are  considered  in  section  2.2.  Therefore 
an  analysis  is  made  on  that  equation. 

Eliminating  pressure,  the  combined  momentum  equation 
becomes 


*2 
d  u 

dye)' 


-\2 

O  V  ,  o  u 

dxc)T  U  dydx 


-  u 


afv 

Sk2 


v 


dy2 


v 


^2 
O  V 

dxdy 


* 


The  result  is  the  same  as  for  the  combined  equations. 


Li 


by' 


bybx 4 


bx' 


dxdy' 


bx 


Replacing  the  derivatives  by  finite  difference  quotients , 
then  substituting  the  approximate  magnitude  of  v,  the 
above  equation  may  be  reduced  to 

^2U  .  d2U  _  d3u  d3u  _  cW_ 

Sy^-r  S^Sx  ay3  aySx2  -  sx  • 

Since  x  is  considered  to  be  infinitely  large,  it  fol¬ 
lows  from  the  continuity  equation  that 


du  _ 
dx 


and 


3v 


/ 


consequently  this  permits  the  combined  momentum  equation  to 
be  reduced  to 

c)2u  =  b^u  _  b<$> 

3y3r  ^y3  3x 


Now  consider  energy  equation  2.1-7.  Neglecting  dissi¬ 
pation  a  priori ,  the  equation  becomes 


c)<t> 

37 


+ 


u 


3<t>  .  d4> 

5?  +  v  5? 


l  /52<t> 

0  5x2  Sy2 


*  Also  justified  on  the  grounds  of  a  constant  surface  tempera¬ 
ture  gradient. 


. 


Again  replacing  the  derivatives  by  finite  difference  quo¬ 
tients  ,  then  applying  the  above  approximate  magnitudes  on 
v  and  x,  it  is  evident  that 


13 


cF 


+ 


1 

a 


Assuming  that  the  remaining  advection  term  is  negligible  the 

* 

boundary  layer  equations  become 


>2 
o  u 

c3y<F 


dy3 


£<t> 


d<t> 

S7 


2.1-9 


2.1-4  Possible  Similarity  Solution 

It  will  now  be  shown  that  a  similarity  solution  exists 
(i.e.  equations  2.1-9  may  be  reduced  to  a  set  of  linear  ordin¬ 
ary  differential  equations)  for  a  large  horizontal  plate 
whose  surface  temperature  gradient  is  a  function  of  time  only. 

The  continuity  equation  implies  the  existence  of  a  dim¬ 
ensionless  stream  function  ijj  —  Tf/(x, y ,t)  such  that 


*  Justified  a  posteriori  (see  Appendix  A) . 


sldxpiJ 


1 


14 


u 


M 

by 


and  v 


bTj/ 

bx 


Introducing  a  similarity  variable 


T)(y/T)  =  ye  (t) 


and  dependent  variables 


^ (y,T)  -  £  (t)  f(t]) 


<t>(x,y,T)  =  C(x,t)  $(tj) 


in  equations  2.1-9,  we  seek  the  unknown  functions  e(r), 

| (t)  and  £ (x,t)  such  that  the  resulting  equations  transform 
to  a  system  of  ordinary  differential  equations  in  F  and  <2>. 
Substituting,  gives 


0  0  I 

-  a1rj  F 


2a^  + 


I  8 

F 


a3$  =  0  2.1-10 


and  $  -  c^r)  $  -  a4  $  -  0 


2 . 1-11 


where 


2 . 1-12 


2.1-13 


2.1-14 


5 
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2.1-15 


Thus  a  similarity  solution  will  exist  provided  it  can  be 
shown  that  all  the  coefficients  a  are  constant. 

It  is  obvious  that  there  are  more  conditions  available 
than  unknowns,  and  therefore  not  all  a's  will  be  arbitrary. 
Proceeding  then  to  the  permissible  forms  of  the  unknowns, 
an  elimination  of  €  between  equations  2.1-12  and  2.1-13 
results  in  the  following  ordinary  differential  equation  which 
may  be  solved  for 
That  is. 


where 


2.1-16 


The  permissible  solutions  are 


y  y  —  1 

i  =  (V  +  V  1  ;  if  - -  *  1  2.1-17a 


7- 


*3  e 


V 


;  if 


?!  ~  1 

7, 


1  ,  2 . l-17b 


noting  that  7^,  A2,  A3,  and  A4  are  arbitrary  constants.  It 
is  important  to  note  that  2.1-17b  is  not  a  separate  solution. 
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but  merely  the  asymptote  of  2.1-17a  as  7  — co  . 


Substituting  back  into  2.1-13  for  |  gives  the  permissible 
form  of  e,  namely 


If  £  is  taken  in  the  form  of  a  product  solution 

C(x/T)  =  M(x)  N(t)  , 

an  inspection  of  equations  2.1-14  and  2.1-15  reveals  that 

M(x)  oC  x  . 

Eliminating  €  between  2.1-14  and  2.1-15,  it  can  easily  be 
shown  that 


N  (t  ) 


*A5t  + 


N(t)  =  A?e 


noting  again  that  A^/  A^,  A^  and  Ag  are  arbitrary  constants. 

In  summary,  a  similarity  solution  does  exist  when  the 
plate  is  very  long  and  inertia  and  advection  terms  may  be 
neglected.  Knowing  the  permissible  forms  of  the  unknown 

functions  €,  |  and  £  and  replacing  7  by  n,  it  can  be  shown 

z 

that 


' 


.  :  i  j>  L.£  rus  »-x : 


■*•  v.:  ■  5  on  >  u. 


If/ 


16t 


n+2 


F  (rj)  , 


and 


<\>  = 


n  , 
XT  $ 
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(n>  0) 


Hence  the  boundary  layer  equations  becomes 


IV 

a  F  +  2r]  F 


I  I  0 


I  I 


4  (n  +  1)  F  -  $  =  0 


2.1-18 


I  I 


0  +2t)$  -  4n  $  -  0  ♦ 


2.1-19 


2.1-18  may  be  integrated  once  immediately.  Satisfying 
the  boundary  condition  on  pressure  gives 


oo 


0  0  0 


I  0 


a  F  +  2r)  F  -  2  (2n  +  3)  F  +  j  <P  dr)  =  0.  2.1-20 

The  remaining  boundary  conditions  2.1-8  transformed  to  the 
T)  domain  ares 


F  (0)  =  F  (0)  =*  $  (0)  -1  —  0 


F  (oo )  =  <$>  (oo )  =  0 


2.2  SOLUTIONS  TO  THE  BOUNDARY  LAYER  EQUATIONS 


2.1-21 


Solutions  of  2.1-19  and  2.1-20  consistent  with  2.1-21 


may  be  found  by  initially  solving  for  the  temperature 


I  :  « 
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distribution,  then  substituting  and  solving  for  the  velocity 
distribution.  Various  mathematical  techniques  are  available 
to  determine  these  solutions?  the  method  of  Frobenius  is 
used  below  and  Laplace  Transforms  were  utilized  in  Appendix 
B. 

2.2-1  Series  Solution  For  Temperature 

Consider  a  solution  of  the  energy  equation 


i  i  i 


2.2-1 


0  +  2t)  0  -  4n  0  =  0 


with  boundary  conditions  0(0)  -  1,  0(oo)  =  0.  Assuming 

this  solution  may  take  the  form 


oo 


o(ri)  =  y 


m+cjD 


m=0 


a  recursion  relation  may  be  found  by  substituting  and 
grouping  coefficients  of  like  powers  of  r|,  namely 


a  .  2  (go  +  r  -  2)  -  4n 

r-2 


a 


(co+r)  (cjd  +  r  -  1) 


r 


The  corresponding  indicial  equation  is  given  by 


a  (D  (od  -  1) 


o 


0. 


19 


Considering  the  zero  root,  it  is  evident  that  a  solution  to 
equation  2,2-1  is 


$(tj)  =  a 


,  .  ~  2  2n  .4 

1  +  2nT)  -  — r-  (1  -  n)  T) 


+  (1  -  n)  (2  -  n)  T]6  .  .  „  . 


2,2-2 


+  a. 


-  i  (1  -  2n)  n3  +  Jo  (1  -  2n)  (3  -  2n) 


-  gig  (1  -  2n)  (3  -  2n)  (5  -  2n)  r )7  . 


The  series  solution  corresponding  to  the  root  a)  -  1 
is  proportional  to  the  second  of  the  two  infinite  series 
obtained  above.  Consequently  the  solution  for  go  -  0  which 
contains  two  arbitrary  constants  is  taken  as  the  complete 
primitive. 

Satisfying  the  boundary  condition  0(0)  =  1  requires 

that  a  =  1.  To  obtain  the  sums  of  the  individual  series 
o 

for  large  values  of  rj,  Euler's  transformation  may  be  used 
to  improve  the  rate  of  convergence.  Utilizing  this  trans¬ 
formation,  the  boundary  condition  at  infinity  requires 
that  (15) 


t. 

20 


2T  (n  +  1) 

“  r(n  +  1/2) 


A  re-arrangement  of  the  series  in  terms  of  Weber  functions 
(15)  permits  the  solution  to  be  given  as 


0(t})  =  22n  r  (n  +  1)  i2nerfc  T), 


2.2-3 


where 


2n  =  0,1,2, 


Plots  for  several  values  of  n  are  illustrated  on  figure  2.2. 


2.2-2  Series  Solution  For  Pressure 


Consider  a  solution  of  the  momentum  equation 


§  -  * 


2.2-4 


with  boundary  condition  p(oo)  =  0.  Writing  the  pressure 
departure  as 

p(x,y,T)  =  XTn  Q^), 


equation  2.2-4  may  be  written  as 

^ 


2.2-5 


This  result  may  also  be  found  from  the  method  of  Laplace 
transforms  as  outlined  in  Appendix  B. 


' 
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FIG.  2.2  TEMPERATURE  PROFILES 
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consistent  with  Q(oo)  =  0.  Substituting  the  solution  for 
temperature,  integrating  and  satisfying  the  boundary  con¬ 
dition  at  infinity  permits  the  solution  to  be  given  as 


$(rj)  =  -22n  T  (n  +  1)  i2n+1  erfc  rj 


2.2-6 


where  2n  =  0,1,2 


The  distribution  of  pressure  departure  is  illustrated  for 
several  values  of  n  on  figure  2.3. 

2.2-3  Series  Solution  For  Velocity 

A  solution  of  equation  2.1-20  is  not  possible  without 
previous  knowledge  of  the  temperature  distribution  $.  Since 
$  is  now  available,  a  substitution  into  equation  2.1-20  leads 
to 


'  1  '  •  •  1  ?n  2n+l 

a  F  +  2r)  F  -  2(2n  +  3)  F  =  -2  T  (n  +1)  i  erfc  i> 


2.2-7 


By  inspection,  the  particular  integral  may  take  the  form 


F  (r\)  oc  i 

r 


2n  +  4 


erfc  rj 


Substituting  into  2.2-7  results  in 


-22nr(n  +  1) 


.  2n  +  4 
i 


erfc  T] 


(1  -  a) 


where 


2n  «  0,1,2 


*5  r 


‘ 


n  =  y/ 


1/2 


FIG.  2.3  PRESSURE  PROFILES 


In  seeking  a  solution  of  the  homogeneous  part  of  2.2-7,  it 
may  be  assumed  that  the  complementary  function  takes  the 
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form 


Fc(rl>  =  Xdm  ’l™  +  7  • 
m  =0 


(do  *  0) 


Hence  by  substituting  and  grouping  coefficients  of  like 
powers  of  T)  a  recursion  relation  may  again  be  found;  that  is 


-d 


d_  = 


r-2 


2  (r  -  2  +  7)  (r  -  3  +  7)  -  (4n  +  6)  (r  -  2  +  7) 


(r  +  7)  (r  -  1  +  7)  (r  -  2  +  7) 


■j  « 


It  follows  that  the  corresponding  indicial  equation  is 

d07  (7-1)  (7  “  2)  =0  . 

Considering  the  zero  root,  a  solution  is  given  by 


FC(<1) 


d  +  d, 
o  1 


+  (2n  +  3)  3  +  (2n  +  3)  (2n  +  1)  5 


3a 


30a 


2.2-8 


+  (2n  +  3)  (2n  +  1)  (2n  -  1)  7  + 

630a 


Tl2  + 


(1  +  n) 

3a 


4  +  4n(l  +  n) 

Tj  T  « 

90a 


2n(l  -  n2)  n8 

7  I  I  I  •  •  •  • 

630a 


The  solutions  corresponding  to  the  non-zero  roots  are  again 
constant  multiples  of  the  series  obtained  in  2. 2-8/ therefore 


Li  r--L  •  ••  •  "  n  i.  •  ‘  y:  -  -  '  I 

' 


•  ••  ,  :  ••  •  •••'  i‘.  .  ?WC1  .j 


•  >  ii.  o  >  x  c 


) 
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2.2-8  is  taken  as  the 
The  total  solution  of 
Fc(t])/  namely 


complete  complementary  solution. 

2.2-7  is  then  the  sum  of  F  (n)  and 

P 


F(t))  =  dQ  +  d1 


T)  + 


2n  +  3 

3a 


T)3  + 


(2n  +  3) (2n  +  1) 

2 

30a 


ti5  + 


(2n  +  3)  (2n  +  1)  (2n  -  1)  7 

*5  T1 

630a 


2.2-9 


+  d„ 


,2  +  SJL+J^L  n4  +  4S.(1  +  .n)  -6 


3a 


90a 


n  - 


2n(l_Z_nj)_  8 

II  •  *  o  • 

630a 


22n  r(n  +  1)  _.2n  +  4 
(1  -  a)  1 


erfc  i). 


Satisfying  the  first  boundary  condition  F(0)  *=  0  requires 
that 


d  = 
o 


r(n  +  1) 

16r  (n  +  3)  (1  -  a) 


and  similarly  satisfying  F  (0)  =  0  requires 


di  = 


r(n  +  1) 


16r  (n  +  |)  (1  -  a) 


The  remaining  boundary  condition,  F  (co )  =  0,  is  found  to 


give 


^  _  -22nr (n  +  1) 

d2  8r  (n  +  2)  (1  -  a) 


■ 

YImikh  -»  (p  )  ^  I 


+  v  lL±jrl  IL  +  ,n£l  +  ep  e^±_ a£  +  p 

“V:  )£ 


.  V  (I  _;  .  J  .  .  '  J  _  iil 

I  .  .  .  .  + 


.  r _i  X , .  dp  ia-+J’  i-r  +  V  I&4M  .  * 


p 

jC9 


»  +  (I  •  :  }  r 

*  - (T- b 


••:  ' t  '  .  .  •;  0  3  (0)1  n<;  Yt;:,aoT  -irs:  .loo  ja^:Yx  .*  :l  ;.n  .  .  . 


•  ' 


’ -i  x  (  cu)  ?.  n  .icy  '  -  nuoc  :  n  .«•  f: 


by  successive  application  of  Euler's  transformation  on  the 
first  two  infinite  series  of  2.2-9.  A  single  differen¬ 
tiation  of  2.2-9  then  gives  the  velocity  distribution  in 
closed  form  as 


F  (r|) 


(n  +  1) 


(n  + 


!> 


(1  -  c) 


L 


1  +  (2n  +  3) 


( 3 


) 


(2n 


i-  3><2^  +  (  a  '4  +  . . . . 

6  ' 


22n  T(n  +  1) 

87  (n  +  2)  (1  -  c) 


lii+nl  f  n  +  4a(i  ±_a )  f  n  'f  + 

L  v rz ■  3  15  -  J 


,2n 


r(n  +  1) 

(1  -  a) 


2n  +  3 


erfc  t) 


from  which  it  follows  that 


F  (t}) 


,2n 


(n  +  1) 


(1  -  o) 


.  2n  +  3  . 

i  erfc 


.  2n  + 
-  i 


^erfc  ^  -* 


i:: 


2.2-10 


The  velocity  distribution  is  illustrated  for  selected  values 
of  n  on  figure  2.4  for  water  and  on  figure  2.5  for  glycerine. 
On  both  diagrams,  the  temperature  profile  (n  =  1)  is  given 


for  reference. 
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2.3  GENERAL  LINEARIZED  SOLUTIONS 


Equations  2.2-3,  2.2-6  and  2.2-10  represent  solutions 
to  the  problem  when  the  surface  temperature  gradient  varies 
as  a  simple  power  of  time.  Due  to  the  linear  nature  of  the 
boundary  layer  equations,  a  more  general  solution  may  be 
obtained  by  superposition. 

Assume  the  surface  temperature  can  be  expanded  in  a 
Maclaurin  series  of  time;  that  is  taking 


x  y 

s  =0 


where 


T“0 


o 


It  is  immediately  apparent  that  the  temperature  distribution 
becomes 


oo 


d>(n)  =  y  22n  r(n  +  1)  bn  i2n  erfc  rj  .  2.3-1 


n  =0 


Substituting  2.3-1  into  2.2-5,  the  pressure  distribution 

for  the  same  arbitrary  surface  temperature  gradient  is 

given  by  ^ 

Q(t])  =  “  /  2  F(n  +  1)  bn  r 

n-0 


erfc  T). 


•fi-Xi  •  *•' 

, 

' 


' 

■  [  -  £  >  •  >•  '  5  ' 


A  similar  substitution  of  <E>  into  2.1-20  may  be  used  to 
show  that  the  velocity  profile  is 
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F  (tj) 


GO 

22n  r(n  +  1)  bn 
n  =0 


. 2n  +  3  _ 

i  erfc 


1  " 


i2n  +  3erfc  (  a  > 

K  fa  y 


CHAPTER  III 


EXPERIMENTAL  WORK 


3.1  DESIGN  OF  APPARATUS 


In  view  of  the  approximations  made  for  the  linearized 
theory  in  Chapter  II,  experiments  were  necessary  to  estab¬ 
lish  their  validity  over  the  range  of  conditions  considered 
for  similarity. 

A  pictorial  view  of  the  apparatus  used  is  given  in 
figure  3.1  and  a  detail  view  in  figure  3.2.  It  consisted 
of  a  25  x  6  x  12  inch  tank  of  transparent  walls  (Perspex) , 
the  base  of  which  was  in  the  form  of  a  flat  opaque  plastic 
(Teflon).  This  plate,  it  was  felt,  would  permit  the  estab¬ 
lishment  of  temperature  gradients  longitudinally  at  the 
surface  with  negligible  conduction  losses  for  the  duration 
of  a  test. 

Differential  heating  of  the  surface  was  achieved  by 
means  of  thermoelectric  cooling  modules  which  were  controlled 
by  a  series  of  Ohmite  adjustable  resistors.  A  typical  mod¬ 
ule  is  shown  in  figure  3.3.  It  is  rated  at  22.2  watts  for  a 
maximum  current  of  thirty  amperes  when  the  hot  side  is  main¬ 
tained  at  27 °C  and  the  cool  side  is  thermally  insulated. 

An  electrically  insulating  film  on  both  sides  eliminated  the 
possibility  of  a  short  circuit.  In  all,  thirty  six  modules 
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FIG.  3.1  GENERAL  VIEW  OF  APPARATUS 


FIG.  3.2  DETAIL  VIEW  OF  TANK 


COLD  END 
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FIG.  3.3  PLAN  VIEW  OF  A  MODULE 


FIG.  3.4  SCHEMATIC  OF  HEATING  CIRCUIT 


HOT  END 


V 


34 


were  needed.  These  were  arranged  in  sets  of  four  laterally 
and  connected  to  the  resistors  as  shown  on  the  schematic 
diagram  3.4.  A  Hewlett  Packard  D.C.  power  supply  rated  at 
forty  five  amperes  and  eighteen  volts  was  utilized  to  pro¬ 
vide  the  necessary  power.  Due  to  the  apparent  linear  be¬ 
havior  of  the  modules  and  a  fine  degree  of  adjustment  on 
the  resistors,  a  wide  range  of  linear  rates  of  heating 
could  be  obtained  at  the  surface  of  the  teflon  plate. 

To  support  the  modules  and  the  tank,  a  26  x  6  7/8  x 
18/32  inch  aluminum  plate  was  selected  on  the  basis  of  its 
high  thermal  capacity  and  its  low  specific  weight.  Consid¬ 
ering  the  plate  as  a  simply  supported  beam  carrying  a  uni¬ 
form  distributed  load,  a  maximum  centerline  deflection  of 
0.003  inches  had  been  calculated  based  on  anticipated  depths 
of  fluid  (three  inches  of  water  or  two  inches  of  glycerine) . 
Additional  stiffness  was  introduced  by  the  four  walls  of  the 
tank  and  further  stiffening  ribs  were  mounted  on  the  outside 
of  the  tank  and  therefore  it  was  felt  that  the  deflection 
of  the  horizontal  surface  was  negligible. 

The  levelling  of  the  horizontal  surface  was  carried  out 
in  two  steps.  First,  a  course  adjustment  was  obtained  with 
screws  on  the  apparatus  stand  and  second,  a  more  precise 
adjustment  was  achieved  with  screws  mounted  directly  on  the 
tank  base  (that  is  the  aluminum  plate) .  This  arrangement 
permitted  an  accuracy  of  +  0.25  minutes  in  the  setting  of 
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the  level  longitudinally  with  the  levelling  device  of  figure 
3,5.  Because  of  its  length,  this  device  could  not  be  placed 
fully  in  a  lateral  position,  hence  an  accurate  setting  of 
the  lateral  level  of  the  teflon  surface  could  not  be  achieved. 
In  view  of  the  isothermal  condition  laterally  and  tempera¬ 
ture  gradients  existing  only  longitudinally,  it  would  seem 
that  this  is  not  critical  (see  Appendix  C) . 


FIG.  3.5  LEVELLING  DEVICE 
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3.1-1  Assembly  of  Plate 

One  of  the  critical  points  in  the  design  of  the  appara¬ 
tus  was  the  flatness  of  the  horizontal  surface.  Initially 
the  supporting  aluminum  plate  was  machined  on  both  sides  to 
a  flatness  of  +  0.001  inches.  Then  the  thermoelectric  cool¬ 
ing  modules  were  bonded  to  the  plate  with  an  Epoxy  resin 
cement.  Following  this,  a  coat  of  teflon  bonding  cement 
was  applied  on  the  opposite  side  of  the  modules,  then  a 
0.070  inch  thick  teflon  sheet  was  clamped  to  the  modules. 

As  a  final  step,  this  "sandwich"  was  machined  only  on  the 
teflon  side  until  the  upper  surface  and  lower  surface  were 
parallel  to  within  +  0.001  inches.  Local  roughness  on  the 
teflon  surface  was  minimized  by  polishing  with  emery  paper 
to  a  smooth  finish. 

3 . 2  INS  TRUMENTATION 


From  the  similarity  analysis  of  Chapter  II,  it  was  pre¬ 
dicted  that  solutions  to  the  governing  equations  are  possible 
provided  a  linear  temperature  distribution  is  imposed  at 
the  surface.  Experimentally,  this  gradient  could  be  obtained 
by  simultaneous  manipulation  of  the  adjustable  resistors  in 
parallel  with  the  modules.  To  determine  the  actual  tempera¬ 
ture  of  the  surface,  nine  evenly  spaced  thermocouples  were 
imbedded  in  the  teflon.  0.020  x  0.025  inch  grooves  were  cut 
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laterally  into  the  surface  using  a  scriber  such  that  the 
thermocouples  could  be  placed  along  isothermals  (noting 
that  the  apparatus  had  been  constructed  to  give  nearly 
isothermal  conditions  traversely) .  Honeywell  Type-J  iron- 
constantan  (0.010  inches  in  diameter)  thermocouples  joined 
by  a  lap  weld  were  used. 

All  of  the  above  mentioned  thermocouples  were  con¬ 
nected  to  a  switching  unit  whose  single  output  was  then 
recorded  on  a  Varian  single  channel  strip  chart  recorder. 

An  actual  speed  of  sixteen  inches  of  paper  per  minute  with 
a  full  scale  deflection  of  one  millivolt  in  0.5  seconds 
permitted  temperature  transients  to  be  recorded  with  a  re¬ 
sponse  of  60°F/sec.;  uncertainty  in  reading  the  strip  chart 
recordings  was  +  0.10°F. 

Temperatures  in  the  water  near  the  Teflon  surface  (that 
is  in  the  boundary  layer)  were  also  measured  by  thermo¬ 
couples.  These  were  identical  to  those  mentioned  previously 
except  a  butt  weld  was  used  to  join  the  two  dissimilar  met¬ 
als.  In  all,  four  thermocouples  were  mounted  in  a  staggered 
manner  on  a  common  stem,  the  position  of  which  could  be 
accurately  controlled  with  a  micrometer  screw  mounted  on 
top  of  the  tank.  A  resolution  of  +  0.0005  inches  in  locat¬ 
ing  thermocouples  from  the  surface  permitted  a  sufficient 
number  of  temperature  recordings  to  be  taken  over  a  rela¬ 
tively  small  distance  (that  is  over  anticipated  boundary 
layer  thicknesses  of  0.10  inches).  These  thermocouple 
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responses  were  recorded  on  two  single  channel  strip  chart 
recorders,  identical  to  that  mentioned  above.  Bulk  fluid 
temperatures  were  determined  by  four  mercury-in-glass  ther¬ 
mometers  which  were  located  at  the  inside  corners  of  the 
tank  such  that  they  would  not  disturb  the  main  flow. 

3.3  EXPERIMENTAL  METHOD 


As  evident  from  the  theoretical  predictions  of  Chapter 
II,  the  linearized  theory  is  applicable  to  fluids  in  a 
large  range  of  Prandtl  numbers.  Since  it  is  neither  desir¬ 
able  nor  possible  to  investigate  all  fluids,  experiments 
were  performed  only  in  water,  glycerine  and  air. 

The  theoretical  predictions  also  embrace  a  wide  range 
of  temporal  variations  of  the  surface  temperature.  One 
which  is  relatively  easy  to  set  up  in  a  laboratory  is  the 
linear  variation,  hence  it  was  decided  to  carry  out  experi¬ 
ments  imposing  linear  surface  temperatures  with  respect  to 
both  time  and  position  along  the  plate.  This  then  fixes  the 
temperature  index  n  at  unity  resulting  in  c  being  the  only 
parameter  for  all  tests. 

3.3-1  Testing  Program 

In  order  to  obtain  a  sufficient  amount  of  experimental 
data,  the  experimental  schedule  was  carried  out  in  two 
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phases.  First,  a  number  of  investigations  were  made  in 
different  fluids  to  determine  the  validity  of  the  approxi¬ 
mations  for  the  linearized  theory.  Then  tests  were  under¬ 
taken  to  determine  the  uncertainty  in  the  main  test  results. 

Main  Schedule: 

(a)  Initially,  experiments  were  made  in  water  for  three 
different  values  of  c  which  not  only  provided  results 
but  also  gave  an  indication  of  the  limits  of  the 
apparatus  and  the  time  taken  for  leading  edge  effects 
to  propagate  along  the  plate. 

(b)  Secondly,  a  test  in  glycerine  was  performed  to  gain 
knowledge  about  the  behavior  of  the  flow  for  at  least 
one  large  Prandtl  number. 

(c)  Finally,  an  investigation  in  air  was  made  to  determine 
whether  or  not  a  radiation  heat  loss  was  appreciable. 

Supplementary  Tests  Schedule: 

(a)  Since  a  change  in  angle  (with  respect  to  the  horizontal) 
of  the  working  surface  may  have  an  effect  on  the  main 
flow,  a  test  in  water  at  c  =  0.34°F/ft.  sec.  was  made 
for  various  angles. 

(b)  A  final  auxiliary  test  was  carried  out  in  water  for 
various  depths  of  fluid  to  determine  the  effect  of 


back  flow. 
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3.3- 2  Preliminary  Tests 

Difficulties  (due  to  the  delicacy  of  the  thermoelectric 
modules)  during  the  initial  trial  runs  led  to  malfunctions 
in  a  few  modules.  Consequently,  the  working  section  had 
to  be  reduced  from  an  initial  length  of  twenty  five  inches 
to  a  final  length  of  fifteen  inches  by  inserting  a  baffle 
at  one  end  of  the  tank.  This  was  found  adequate  for  test 
purposes . 

At  an  advanced  stage  in  the  experiments  with  water, 
the  surface  thermocouples  gave  erratic  readings.  Initially, 
those  at  the  hot  end  of  the  plate  became  unreliable  and  ev¬ 
entually  all  of  the  surface  thermocouple  readings  became 
irregular.  This,  it  is  believed,  can  be  attributed  to 
differential  expansion  between  the  thermocouples,  cement 
and  teflon  ultimately  leading  to  the  separation  of  the 
junctions.  The  probe  thermocouples  were  therefore  used  in 
both  setting  boundary  conditions  and  in  the  measurement  of 
temperatures  in  the  boundary  layer. 

3.3- 3  Test  Procedure 


Considering  a  test  in  water  the  procedure  (which  was 
the  same  for  all  runs)  consisted  of  the  following  steps. 
Utilizing  one  of  the  probe  thermocouples ,  the  surface  temp¬ 
erature  was  adjusted  (by  varying  the  power  supply  to  sets 
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of  modules)  such  that  it  varied  both  linearly  with  respect 
to  position  along  the  plate  and  time.  Typical  surface  temp¬ 
erature  distributions  at  various  times  are  illustrated  on 
figure  3.6  for  c  =  0.49°P/ft.  sec.  The  apparatus  was  then 
allowed  to  stand  undisturbed  at  room  temperature  until  it 
was  considered  to  be  isothermal  throughout.  The  probe  thermo¬ 
couple  was  then  located  in  the  boundary  layer  region  and  a 
strip  chart  recorder  was  started.  Switching  on  the  power 
supply  immediately  after  this,  a  "blip"  appeared  on  the  strip 
chart  providing  an  excellent  timing  mark.  The  temperature 
was  then  continuously  recorded  until  a  definite  change  in  the 
linear  behavior  was  observed  (indicating  the  onset  of  a  lead¬ 
ing  edge  effect) . 

Since  each  temperature  was  measured  at  only  one  point 
at  a  time  all  of  the  curves  were  superimposed  until  the 
"blips"  coincided  hence  giving  reference  to  the  same  temp¬ 
oral  origin.  The  reproducibility  of  the  curves  was  found 
to  be  good. 
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FIG.  3.6  TYPICAL  SURFACE  TEMPERATURE  PROFILE 


CHAPTER  IV 


EXPERIMENTAL  RESULTS 


The  results  of  tests  performed  in  water,  glycerine 
and  air  are  presented  in  this  chapter. 

A  spot  check  obtained  by  immersing  the  thermocouple 
in  firstly  a  mixture  of  ice  and  water  and  then  boiling  water 
the  readings  of  which  were  taken  on  a  potentiometer,  read 
the  same  temperature  as  an  immersed  mercury-in-glass  thermo¬ 
meter.  Hence,  assuming  a  linear  correspondence  between  the 
strip  chart  recorder  reading  and  the  temperature  at  the 
thermocouple  junction,  all  intermediate  changes  in  milli¬ 
volt  output  were  converted  to  changes  in  temperature  by 
using  the  manufacturer's  calibration  chart. 

4 . 1  TESTS  IN  WATER 

Having  used  ordinary  tap  water  with  no  special  prep¬ 
aration  other  than  natural  de-aeration  over  a  twenty  four 
hour  period,  comparisons  of  experimental  results  obtained 
with  theoretical  temperature  distributions  are  given  in 
figures  4.1,  4.2  and  4.3.  A  sample  calculation  of  the  con¬ 
version  of  the  strip  chart  recordings  may  be  found  in  Appen¬ 
dix  D.  Figure  4.1  illustrates  typical  variations  of  tempera 
ture  with  time  at  various  positions  from  the  surface  (for 
c  =  0.34 °F/ ft.  sec). 
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FIG.  4.1  VARIATION  OF  TEMPERATURE  WITH  TIME 


In  view  of  the  very  small  temperature  increases  as¬ 
sociated  with  the  edge  of  the  boundary  layer  where  experi¬ 
mental  errors  might  be  very  large,  no  results  are  given 
for  this  range.  An  estimated  thickness  of  the  thermal 
boundary  layer  is  0.12  inches.  Due  to  inaccuracies  in  esti 
mating  small  distances  from  the  surface,  temperature  meas¬ 
urements  were  not  made  within  0.010  inches  from  the  surface 

Agreement  in  figure  4.1  is  mostly  good  for  small  times 
but  for  large  times ,  when  the  effect  of  inertia  is  expected 
to  become  more  important,  the  discrepancy  increases.  The 
agreement  is  poorest  near  y  *  0.040  inches,  in  the  region 
of  velocity  maximum. 

Figure  4.2  illustrates  the  variation  of  temperature 
with  height  at  various  times  for  the  same  conditions  as  in 
figure  4.1.  Near  the  surface  and  at  the  edge  of  the  bound¬ 
ary  layer  where  viscous  forces  tend  to  dominate  agreement 
is  good.  The  effect  of  inertia,  if  present,  would  assume 
its  greatest  significance  in  the  region  of  maximum  velocity 
Due  to  inertia  the  flow  might  be  expected  to  be  retarded 
hence  experimental  temperatures  would  lie  above  theory. 

This  certainly  seems  to  be  evident  in  the  maximum  velocity 
region  of  figure  4.2  where  the  discrepancy,  although  not 
large,  is  greatest.  Temperature  recordings  are  not  shown 
for  larger  times  than  sixteen  seconds  due  to  an  onset  of  a 
leading  edge  effect  as  discussed  in  Chapter  5,  section  5.2. 

The  non-dimensional  temperature  profiles  for  two  tests 
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HEIGHT  (INS) 


FIG.  4.2  VARIATION  OF  TEMPERATURE  WITH  HEIGHT 
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are  compared  with  the  linearized  theory  (n  =  1)  in  figure 
4.3.  A  third  test  (c  =  0.49°F/ft.  sec.)  gave  erratic  date 
due,  it  is  thought,  to  incipient  instability:  it  is  not 
included  in  the  results  but  is  discussed  in  section  5.4. 

The  velocity  curve  (n  =  1)  is  given  for  reference  purposes. 
The  points  which  were  shown  in  figures  4.1  and  4.2  corres¬ 
ponding  to  c  =  0.34°F/ft.  sec.  again  indicate  fair  overall 
agreement.  Possible  experimental  errors  near  the  surface 
in  $  and  rj  are  +  0.04  and  +  0.02  respectively,  due  mainly 
to  inaccuracy  in  measuring  small  distances  from  the  surface. 
In  the  maximum  velocity  region,  experimental  errors  in  $  and 
T)  are  about  +  0.04  and  +  0.02  respectively,  indicating  a 
possible  effect  of  inertia  although  not  confirming  it  defini¬ 
tely.  Possible  errors  in  $  of  +  0.02  and  in  r)  of  +  0.02 
near  the  edge  of  the  boundary  layer  are  mainly  due  to  in¬ 
accuracy  in  detecting  small  temperature  increases  associated 
with  these  points.  The  results  for  c  =  0.27°F/ft.  sec.  show 
good  agreement  over  the  entire  profile  with  the  exception  of 
the  largest  value  of  tj.  A  possible  experimental  error  in 
of  +  0.02  in  this  region  again  is  a  result  of  inaccuracy  in 
estimating  small  temperature  increases. 

A  fourth  test  (c  =  0.34°F/ft.  sec.)  performed  at  a 
distance  of  seven  inches  from  the  leading  edge  was  in  simi¬ 
lar  agreement  with  the  theory.  Since  no  additional  informa¬ 
tion  about  the  effect  of  inertia  could  be  obtained,  it  was 


not  included  in  this  thesis. 
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FIG.  4.3  COMPARISON  OF  TEMPERATURE  PROFILES  (WATER) 
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4.2  TEST  IN  GLYCERINE 


The  non-dimensional  temperature  profile  for  one  test 
in  glycerine  is  compared  with  the  linearized  theory  in 
figure  4.4.  For  large  Prandtl  numbers,  where  the  thermal 
boundary  layer  is  much  thinner  than  the  momentum  boundary 
layer,  it  might  be  expected  that  the  viscous  forces  domin¬ 
ate  the  flow.  In  view  of  the  overall  agreement  it  certain¬ 
ly  seems  to  be  true  for  this  test  (c  =  0.49°F/ft.  sec.), 
hence  suggesting  that  the  linearized  theory  is  also  adequ¬ 
ate  for  predicting  temperature  profiles  in  glycerine  during 
leading  edge  free  conditions.  Again,  possible  experimental 
errors  in  $  and  T)  are  +  0.030  and  +  0.015  respectively  near 
the  surface  (small  rj)  ,  whereas  for  larger  values  of  T}  the 
errors  in  4>  and  t)  become  +  0.20  and  +  0.020  respectively. 
Discrepancies,  although  not  large,  seem  to  be  greatest  near 
the  surface  which  can  be  attributed  to  the  error  in  estimat¬ 
ing  small  distances. 

The  velocity  profile  (n  =  1)  is  also  given  to  reveal 
the  monotonic  nature  of  the  curve  within  the  thermal  bound¬ 
ary  layer. 

4.3  TEST  IN  AIR 


It  has  been  shown  (16)  that  the  effects  of  radiation 
losses  in  air  for  free  convection  on  an  inclined  plate  are 
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FIG.  4.4  COMPARISON  OF  TEMPERATURE  PROFILES  (GLYCERINE) 
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quite  high.  To  determine  a  quantitative  indication  of  this 
effect  for  the  present  problem,  a  test  was  carried  out  in 
air  for  c  =  0.52°F/ft.  sec.  The  experimental  results  are 
compared  with  the  linearized  theory  on  the  non-dimensional 
plot  of  figure  4.5.  Since  overal  agreement  was  good  for 
both  the  water  and  glycerine  tests  and  because  the  linear¬ 
ized  theory  has  neglected  radiation  losses,  it  would  seem 
reasonable  that  the  disagreement  might  be  due  to  this  effect. 
Possible  experimental  errors  in  0  and  r\  of  +  0.030  and  +  0.020 
respectively  and  an  estimated  response  time  of  0.5  seconds 
for  the  thermocouple  junction  in  air  would  not  account  for 
the  disagreement.  If  radiation  losses  are  present,  all 

thermocouple  readings  would  be  expected  to  be  high  and  hence 

★ 

the  extrapolation  technique  would  predict  a  higher  surface 
temperature  (at  a  given  time)  than  actually  existed  on  the 
plate.  In  plotting  figure  4.5,  all  temperatures  were  nor¬ 
malized  with  this  surface  temperature  thereby  suggesting 
that  experimental  results  would  be  below  the  theory. 

On  the  basis  of  this  test  radiation  losses  appear  to 
be  about  50  per  cent  of  the  total  heat  loss. 
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see  Appendix  D. 
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FIG.  4.5  COMPARISON  OF  TEMPERATURE  PROFILES  (AIR) 
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CHAPTER  V 


DISCUSSION 


Chapter  IV  has  presented  the  experimental  results 
of  the  main  tests  with  a  discussion  of  errors.  This 
chapter  considers  the  validity  in  the  main  test  results. 

5 . 1  LATERAL  EDGE  EFFECTS 


A  tacit  assumption  in  both  the  theoretical  and  experi¬ 
mental  work  is  that  the  lateral  edges  of  the  heated  hori¬ 
zontal  surface  are  sufficiently  far  removed.  These  edges 
promote  a  boundary  layer,  the  growth  of  which  could  inter¬ 
fere  with  the  main  flow  and  destroy  the  two-dimensionality 
of  the  problem. 

Although  the  intensity  of  the  effects  produced  by  the 
side  edges  will  vary  for  different  fluids,  they  will  in  gen¬ 
eral  be  minimized  if  the  edges  are  designed  for  the  worst 
conditions.  To  eliminate  any  direct  buoyancy  effects,  the 
edges  were  maintained  at  a  nearly  isothermal  condition. 
Anticipated  thicknesses  of  the  main  boundary  layer  in  air 
(0.30  inches)  provided  a  design  criterion  for  locating  the 
side  walls.  If  the  thickness  of  the  side  wall  boundary 
layer  is  taken  to  be  of  the  same  order  of  magnitude  as 
the  thickness  of  the  main  flow,  six  inches  between  the  walls 
would  seem  sufficient  to  remove  side  wall  effects. 
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Experimentally,  these  effects  may  be  revealed  by  a 
temperature  traverse  laterally  at  a  fixed  distance  from 
the  leading  edge  and  above  the  horizontal  surface.  Such 
a  traverse  was  made  in  both  water  and  air  with  the  probe 
thermocouple  located  at  the  furthest  possible  position 
from  the  leading  edge  (12  inches) .  The  results  of  the 
experiment  in  water  are  given  in  figure  5.1,  where 
c  =  0.34 °F/ ft.  sec.,  for  various  times.  It  is  noticed 
that  a  boundary  layer  formed  near  the  edges  but  it  did 
not  seem  to  affect  the  central  region  (about  4  inches  in 
width) .  During  the  preliminary  tests ,  the  tank  was  found 
to  be  inadequately  sealed  at  the  base,  hence  it  was  nec¬ 
essary  to  spread  silicone  rubber  cement  around  the  edges. 
As  a  result,  this  would  account  for  the  large  region  of 
influence  at  the  side  walls.  The  test  in  air  revealed  a 
similar  temperature  traverse. 


5.2  LEADING  EDGE  EFFECTS 


Linearization  of  the  flow  equations  was  possible  by 
considering  a  horizontal  surface  of  infinite  extent,  a  con¬ 
dition  which  completely  neglects  leading  edge  effects.  It 
is  therefore  doubtful  if  the  linearized  theory  applies  at 
a  point  in  the  system  where  leading  edge  effects  occur  for 
at  such  a  point  the  effects  of  inertia  and  advection  may  be 
significant.  Simulation  of  the  infinite  condition  would 
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FIG.  5.1  LATERAL  TEMPERATURE  TRAVERSE 
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entail  the  use  of  an  impractical  apparatus,  hence  a  finite 
plate  (14  inches  long)  was  used  in  combination  with  a  de¬ 
tection  method  for  predicting  the  onset  of  a  leading  edge 
effect.  Preliminary  runs  revealed  a  distinct  change  in 
the  behavior  of  the  thermocouple  output  from  the  strip  chart 
record  of  temperature  versus  time.  Consequently,  the  time 
taken  before  this  change  was  adopted  for  predicting  the 
onset  of  leading  edge  effects. 

The  lengths  of  test  periods  differ  considerably  in 
water  (about  20  seconds) ,  air  (about  12  seconds)  and  gly¬ 
cerine  (about  60  seconds) :  only  the  tests  in  water  are 
discussed  here.  By  choosing  a  station  about  two-thirds 
of  the  test  length  downstream  of  the  leading  edge  it  was 
found  possible  to  obtain  a  sufficiently  long  test  period 
and  yet  avoid  any  possible  disturbance  caused  by  backwash 
from  the  trailing  edge.  A  survey  of  the  time  taken  before 
the  change  in  the  thermocouple  output  was  then  made  for 
various  heights  above  the  heated  surface.  This  is  shown 
in  Table  I  for  three  tests.  As  would  be  expected,  it  is 

DISTANCE  PERIOD  OF 

FROM  "ONE -DIMENSIONAL"  BEHAVIOR 

SURFACE  (SECS) 

(INS)  c  =  0.27  c  =  0.34  c  =  0.49 


0.010 

25 

21 

20 

0.020 

23 

18 

16 

0.035 

23 

18 

15 

0.055 

23 

18 

16 

0.095 

24 

20 

19 

TABLE  I  ONSET  OF  LEADING  EDGE  EFFECT 
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evident  that  the  time  decreases  as  the  rate  of  change  of 
surface  temperature  gradient  (c)  increases.  Also  notice¬ 
able  is  a  definite  minimum  with  respect  to  height  for  each 
series.  This  is  consistent  with  the  shape  of  the  velocity 
profile  which  might  suggest  a  minimum  leading  edge  free 
period  near  the  velocity  maximum.  In  view  of  the  absence 
of  a  sharply  pronounced  minimum  it  appears  that  the  onset 
of  the  leading  edge  effect  could  be  estimated  quite  accur¬ 
ately  from  the  bulk  velocity. 


5.3  DIRECT  EFFECT  OF  GRAVITY 


For  laminar  natural  convection  flow  on  a  vertical 
plate  the  generating  body  force  acts  parallel  to  the  plate 
whereas  on  a  horizontal  surface  it  acts  normal  to  plate. 

In  the  latter  case  the  body  force  still  generates  the 
flow  but  only  indirectly.  On  an  inclined  plate,  compon¬ 
ents  acting  both  parallel  and  normal  to  the  plate  would 
each  govern  the  flow  behavior  but  for  small  angles  from 
the  horizontal  little  effect  due  to  the  former  might  be 
expected. 

A  series  of  tests  were  made  in  water  by  holding  c 
constant  and  then  recording  temperature  profiles  for  various 
plate  angles.  The  results  are  given  in  figure  5.2a  for 
downward  inclinations  and  in  figure  5.2b  for  upward  inclina¬ 
tions,  where  both  inclinations  are  relative  to  the  leading 
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0  8  16  24 


TIME  (SECS) 

(a)  INCLINED  DOWN 


TIME  (SECS) 

(b)  INCLINED  UP 

FIG.  5.2  VARIATION  OF  TEMPERATURE  WITH  INCLINATION 


■ 


59 


edge.  In  view  of  possible  experimental  errors  in  record¬ 
ing  temperatures  and  time  and  an  uncertainty  of  +  0.25 
minutes  in  inclination,  it  may  be  concluded  that  no  sig¬ 
nificant  change  was  observed  within  an  angle  of  +  15 
minutes  from  the  horizontal. 

5.4  STABILITY 

Two  layers  of  fluid  of  different  densities  with  the 
heavier  above  the  lighter  are  in  a  state  of  metastable 
equilibrium.  Extensive  studies  (3)  have  been  made  on  a 
meteorological  scale  where  it  was  found  that  a  layer  of 
fluid  heated  from  below  would  become  unstable  depending 
on  a  critical  parameter  (the  Rayleigh  number) .  Superposi¬ 
tion  of  a  steady  flow  along  the  horizontal  surface  in  con¬ 
trast  has  the  effect  of  promoting  stability  (17) . 

The  Rayleigh  number  associated  with  thermal  instabi¬ 
lity  in  a  quiescent  layer  of  fluid  for  free  upper  and  lower 
surfaces ,  or  for  a  free  upper  surface  and  a  rigid  lower 
surface  and  for  the  case  of  rigid  upper  and  lower  surfaces 
has  been  defined  as 

_  ftgh3A8 
Ra  "  kv 

where  h  is  the  boundary  layer  thickness  and  A9  is  the 
temperature  difference  between  the  two  surfaces.  Choosing 
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AO  as  the  excess  temperature  near  the  probe  thermocouple, 
the  Rayleigh  number  was  calculated  at  the  end  point  of 
each  test  discussed  in  Chapter  4.  Table  II  summarizes 
the  results.  The  critical  Rayleigh  number  for  a  free  iso¬ 
thermal  upper  surface  and  a  rigid  isothermal  lower  surface 
is  about  1100(18),  hence  Table  II  suggests  that  only  one 
water  test  (c  =  0.49°F/ft.  sec.)  may  have  been  incipiently 
unstable.  This  seems  to  be  in  agreement  with  an  observed 
temperature  oscillation  near  the  end  point  of  that  test. 


Glycerine 

c  =  0.49 

Ra  =  11 

Air 

0.52 

20 

Water 

0.27 

960 

Water 

0.34 

1300 

Water 

0.49 

1700 

TABLE  II 

FINAL  RAYLEIGH  NUMBERS 

MISCELLANEOUS 

OBSERVATIONS 

Observation  on  finite  suspended  particles  in  the  water 
tests  indicated  the  existence  of  a  region  where  viscous 
effects  dominate  near  the  horizontal  surface.  From  these 
observations  it  was  possible  to  estimate  the  bulk  velocity, 
about  2  cm. /sec. ,  for  c  =  0.34°F/ft.  sec.  By  a  similar 
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technique  in  glycerine,  (with  a  dye  streak  in  the  lateral 
direction)  flow  was  also  observed  resulting  in  an  estimated 
bulk  velocity  of(c  =  0.52°F/ft.  sec.)  about  0.25  cm. /sec. 

In  view  of  the  thin  boundary  layer  thicknesses ,  quantita¬ 
tive  measurements  by  these  techniques  were  not  possible. 

Another  aspect  which  might  have  influenced  experimental 
results  appreciably  was  the  formation  of  air  bubbles  in 
water  on  all  surfaces.  These  were  noted  to  be  of  the  order 
of  the  thermocouple  diameter  (0.010  inches).  To  overcome 
this  problem,  each  time  new  tap  water  was  used  it  was 
heated  slightly  and  then  allowed  to  de-aerate  over  a  twenty- 
four  hour  period  hence  eliminating  all  observable  entrained 
air. 

A  final  aspect  which  may  have  influenced  the  experi¬ 
mental  results  is  the  possibility  of  "back  flow"  depending 
on  the  depth  of  fluid.  To  determine  the  effect  of  depth 
of  fluid  on  the  temperature  distribution,  profiles  were  re¬ 
corded  for  various  depths  of  water  with  the  probe  thermo¬ 
couple  located  at  0.030  inches  from  the  surface  and  seven 
inches  from  the  leading  edge  (c  =  0.32°F/ft.  sec.).  Figure 
5.3  illustrates  the  results.  It  is  noticed  that  thermal 
instability  was  not  detected  and  in  view  of  the  insignifi¬ 
cant  disagreement  the  depth  of  fluid  does  not  seem  to  in¬ 
fluence  temperature  profiles. 


, 

;.cL  a&oq  -ox:  s> .  tv  oc/pimiot 

5  3  rtc  :  -  i  -*  q  .  ...  .vi 

•  :cr  .  . 

rf  J  iv? 


62 


TIME  (SECS) 


FIG.  5.3  VARIATION  OF  TEMPERATURE  WITH  TIME 


FOR  VARIOUS  DEPTHS  OF  WATER 
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CHAPTER  VI 


CONCLUSIONS  AND  RECOMMENDATIONS 


6.1  CONCLUSIONS 


This  thesis  has  described  a  theoretical  and  experi¬ 
mental  investigation  of  laminar  natural  convection  in  the 
unsteady  state  on  a  large  horizontal  surface.  The  linear¬ 
ized  theory  developed  for  power  law  variations  (with  respect 
to  time)  of  the  surface  temperature  gradient  was  well  sup¬ 
ported  by  the  experimental  results.  In  view  of  the  range 

of  fluids  tested,  that  is,  those  within  Prandtl  numbers  of 

3 

3.0  and  about  12  x  10  ,  this  indicates  that  the  inertia  and 
advection  terms  are  not  significant  for  non-metallic  liquids 
Although  the  theoretical  predictions  appear  to  be  valid  in 
the  range  of  low  Prandtl  numbers  (liquid  metals) ,  it  remains 
to  be  shown  whether  the  linearized  theory  is  applicable 
there. 

A  comparison  of  experimental  and  theoretical  tempera¬ 
ture  profiles  for  water  revealed  good  overall  agreement. 
Discrepancies,  although  never  large,  appear  to  be  greatest 
near  the  region  of  the  velocity  maximum  and  this  is  consist¬ 
ent  with  the  possible  effect  of  inertia  as  felt  through  the 
velocity  maximum.  A  similar  comparison  for  tests  in  glycer¬ 
ine  also  revealed  good  overall  agreement  with  no  large  devia 
tions  in  any  region  of  the  velocity  profile.  This  is 
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consistent  with  the  expected  dominance  of  viscosity  at 
large  Prandtl  numbers.  An  experiment  performed  in  air  was 
not  in  good  agreement  with  the  linear  theory.  In  view  of 
a  consistent  disagreement,  radiation  losses  seem  to  be 
significant. 

Differences,  between  the  temperature  profiles  ob¬ 
tained  for  various  angles  of  the  plate  from  the  horizontal 
and  the  temperature  profiles  obtained  at  the  horizontal 
position,  were  found  to  be  insignificant  in  the  range  of 
plus  or  minus  fifteen  minutes  of  slope.  Tests  performed 
in  water  under  various  depths  of  fluid  revealed  a  negli¬ 
gible  difference  in  the  temperature  profiles  for  depths  of 
fluid  less  than  three  inches.  Existence  of  flow  was  dem¬ 
onstrated  by  observing  small  particles  suspended  in  water 
and  a  dye  streak  in  glycerine. 

The  results ,  and  hence  the  conclusions ,  are  only  valid 
for  circumstances  in  which  the  effect  of  a  leading  edge  can 
be  ignored,  in  which  the  side  edges  are  sufficiently  far 
removed  and  only  under  moderate  temperate  changes. 

In  view  of  the  consistent  overall  agreement,  the  lin¬ 
earized  theory  would  be  expected  to  accurately  predict  the 
gross  characteristics  of  the  boundary  layer,  particularly 
the  heat  transfer  relations. 
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6.2  APPLICATIONS 

Although  only  the  heated  plate  facing  upwards  has  been 
studied  in  this  thesis,  a  slight  alteration  will  result  in 
similar  expressions  for  a  heated  plate  facing  downwards. 
Cooled  plates  either  facing  up  or  down  may  also  be  analyzed 
on  the  same  basis;  only  the  boundary  conditions  on  the  sur¬ 
face  temperature  gradient  are  altered. 

Applications  of  these  "gravity"  flows  are  most  likely 
realized  on  a  micrometeorological  or  oceanographic  scale; 
in  particular  in  the  study  of  ocean  currents  or  sea  breezes. 
These  problems  are  more  difficult  than  the  present  developed 
linearized  theory  would  suggest,  but  it  might  be  used  as  a 
zeroth  approximation.  In  the  field  of  civil  engineering, 
it  might  be  used  in  the  study  of  density  currents  in  re- 
sevoirs  and  lakes  where  the  mass  transfer  problem  is 
analogous  to  the  heat  transfer  problem. 

6 . 3  RECOMMENDATIONS 


Conclusions  regarding  the  effects  of  inertia  and 
advection  have  been  reported  from  experimental  temperature 
profiles  only.  A  more  appropriate  method  would  be  a  direct 
study  of  the  experimental  velocity  profile.  Complexity 
due  to  small  velocities  in  the  unsteady  state  makes  it 
impossible  at  the  present  time  to  devise  a  direct  velocity 
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measuring  instrument.  It  may  be  feasible  to  use  a  fibre 
anemometer  but  current  steady  state  devices  are  impractical 
for  unsteady  state  conditions.  Indications  of  bulk  velo¬ 
cities  are  attainable  by  observing  suspended  particles  in 
water  and  also  by  timing  the  movement  of  a  dye  streak  in 
glycerine  between  two  points  on  the  horizontal  surface. 

To  justify  the  validity  of  the  approximations  in  the 

linearized  theory  over  the  full  range  of  Prandtl  number, 

* 

the  experimental  work  should  be  extended  to  the  liquid 
metals.  This  would  involve  the  re-designing  of  a  larger 
tank,  since  times  for  the  propagation  of  leading  edge  effects 
are  much  shorter.  Further,  the  solutions  should  be  ex¬ 
tended  to  include  the  non-linear  terms  for  these  might  be¬ 
come  significant  at  low  Prandtl  numbers. 

Laying  thermocouples  beneath  the  surface  was  found  to 
be  impractical  (malfunction  due  to  differential  expansion) 
whereas ,  thermocouples  placed  on  top  of  the  surface  worked 
well.  In  view  of  the  linear  temperature  gradient  near  the 
wall,  the  actual  surface  temperature  was  obtained  by  linear 
extrapolation. 
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APPENDIX  A 


VALIDITY  OF  APPROXIMATIONS 

Consider  the  validity  of  the  approximations  made  in 
Chapter  II  for  a  test  in  water  (c  =  0. 34 °F/ft. sec. ) . 

First  examine  the  orders  of  magnitude  of  the  combined 
momentum  equation 

c^u  _  c^u  d<t>  ,  , 

5y5r  -  ^3  -  5S  •  A_1 

If  the  problem  is  to  be  properly  posed  both  the  unsteady 
and  "driving"  term  must  remain.  It  therefore  remains  to 
be  shown  that  the  viscous  term  is  of  the  same  order  of  mag¬ 
nitude. 

Given  X  =  1.25  ft.,  Y  =  0.0084  ft.,  U  =  0.065  ft/sec., 

X  =  Y  =  0.0095  ft.,  U  =  0.0011  ft/sec., 

o  o  o 

and  t  =  9.0  sec. , 

a  substitution  into  A-l  leads  to 

600/ 1  =  85  -  0.26  ©x  , 

noting  that  the  value  of  U  was  obtained  as  mentioned  in 
Chapter  V,  section  5.5.  Examining  the  relative  orders  of 
magnitude,  it  seems  reasonable  that  each  term  in  A-l  is  of 


the  same  order. 


. 
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Now  consider  approximations  of 
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d<t>  d<t>  _  1  f 

57  +  u  57  -  a  V 


527 

5y2 * 5 *- 


A-2 


Given  a  to  be  about  5.9,  a  substitution  into  A-2  gives 


310/ 1  +  16  =  8 


A©  , 


A-3 


where  A©  has  been  taken  outside  the  bracket  since  it  is  a 
common  factor  when  A-2  is  written  in  terms  of  finite  dif¬ 
ference  quotients.  Examining  A-3  it  appears  that  advection 
may  not  be  neglected  but  since  Chapter  II  considers  an  in¬ 
finite  plate  and  since  x  may  be  any  value  without  changing 
the  other  two  terms,  it  seems  justified  to  neglect  advection. 

Consider  the  validity  of  neglecting  dissipation,  namely 


+  2 


/5u 

dvYl 

Is? 

57/ 

2 

(du/dy)  is  the  only  large  term  in  the  bracket,  therefore 

the  magnitude  of  the  dissipation  term  is 

5  x  lcT7  x  67  . 

As  a  result,  it  seems  justified  to  neglect  dissipation  when 
compared  to  A-3.  It  is  important  to  note  that  near  the  criti¬ 

cal  point  of  a  fluid  where  f3  becomes  large  dissipation  may 
not  be  neglected. 


A 


:  h 


3  - 


s  ai  di  ©on  &  /  irf3  o  >iad  c  '■  j  •’'•• 

-lib  .  ic  i  •  ti-  .c  r  • 

■  >.tdo9vb«  d.i  2.2  a -i  £-A  grin  -  ft  d 

-Hi  ns  sisbigfioD  II  isdqsrfD  sonla  dLrd  bsdoslpsn  sd  don 

. 

;±do©vbs  doslpan  od  b©I3:idaut  bih998  di  ,8fln93  owd  d»rfdo  9rfd 
:9/nsn  ^noidsqiaaib  gnidoalpan  io  Y-iifr-^6'’  ^rfd 


9io  ^i  >rfd  N  -  odd  f  i  •  j..  r  .1  \  '  o  ;  ■ 


mo.  ttol  i  n 


•  3- -  '  >  i  •••  -  d  ■•'•••  g-  1  - . .  >  ’  "  :; 

Yfc  ••:  aoidBQI  taxb  Sr-U  ■  . 


APPENDIX  B 


SOLUTIONS  BY  LAPLACE  TRANSFORMS 


Consider  a  solution  of  the  conduction  equation 


B-l 


with  boundary  conditions  <t>  (o)  =  XTn,  <0(co)  =  0.  The  Laplace 
transform  of  equation  B-l  with  respect  to  the  variable  t  is 


s<0 


1 

a 


/ 


with  boundary  conditions  <0(0,5)  = - ^  ^  ,  4>(oo,s)  =  0 

s 

where  n^>  -1. 

Solving, 

4  =  xr(n  +  D  e-  /as  y 


the  inverse  of  which  is 


<0  =  xtR  22n  f(n  +  1)  i2n  erfc  ly, 


l/2-i 


B-2 


where  2n  =  0 , 1 , 2 ,  . 

Hence,  substituting  B-2  into  the  y-momentum  equation  2.1-9, 
integrating,  and  satisfying  the  condition  on  pressure  at 
infinity  gives 
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P  = 


22n  r(n+l) 


B-3 


To  determine  the  velocity  distribution,  consider  a  solution 
of  the  x-momentum  equation 


c)u  _ 

St  " 


B-4 


consistent  with  u(o)  =  u(oo)  *  0. 

Upon  substitution  of  B-3  into  B-4,  the  Laplace  transform  of 
B-4  with  respect  to  t  using  the  initial  condition 
u(x,y)Ts=0  =  0  yields 


su 


r(n+l)  e“  fe*  y 

1/2  n  +  3/2 
a  '  s  ' 


B-5 


with  u(o,s)  =  u(oo  ,s)  =  0. 

By  inspection  B-5  has  a  particular  integral  given  by 


-  J?  r  (nH-l)  e~  y  . 

up  50  1/2  _n  +  3/2  ' 


s 


substituting  into  B-5  results  in 


S  =  r  (n+l)  el  {p* 
P  sn  +  5/2 


' 


, 
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Seeking  the  complete  solution  in  the  form 


u  =  Ae 


/s~  y  +  Be-/a  y  +  r(n+I) 


-  /as 


a-oa1^2  ' 


B-6 


u  (co  ,s) 
u (o,s) 


0  requires  that  A  *=  0  and 
0  requires 


B  = 


(1-a) a  '  s 


L(P±U _ _  . 

1/2  n  +  5/2 


Hence,  the  inverse  of  B-6  becomes 


u  =  Etei-U  (4T)n+3/2 
(l-a)a1/Z 


i2n+3  erfc  y/(~ 


,1/2 


i2n+3  erfc  y/(4t)V2] 


which  agrees  with  the  expression  2.2-10. 
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APPENDIX  C 


INFLUENCE  OF  LATERAL  VELOCITIES 


If  the  Teflon  plate  were  inclined  laterally,  a  com¬ 
ponent  of  velocity  in  this  direction  would  be  possible. 
The  purpose  of  this  appendix  is  to  determine  the  order  of 
magnitude  of  this  component  hence  permitting  a  comparison 
to  be  made  with  the  velocity  of  the  main  flow. 

Consider  the  inclined  isothermal  plate  of  figure  C.l 
consistent  with  the  boundary  conditions. 


FIG.  C-l  COORDINATE  SYSTEM 


0  (o)  -  Bt  =  U  (o)  =  0 
9  (oo )  =  U(oo)  =  Pd(oo )  =  0. 


Applying  the  same  approximations  as  in  Chapter  II, 
the  governing  boundary  layer  equations  become: 


du  ,  dV  _ 
Sx  +  ” 


0 


f3g©  (cos  x) 


V 
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0  =  f3g0  (sin  X)  "  “  -§y 


c)9 

St 


ic 


Considering  a  solution  of  the  energy  equation ,  the 
Laplace  transform  with  respect  to  the  variable  t  is 


IC 


d20 

dY2 


/ 


with  boundary  conditions  0(o,s)  =  — j  and  0  (00 ,s)  =  0. 

s 

I 

Solving, 


s 


the  inverse  of  which  results  in  the  temperature  distribution 

9  =  B(4t)  i2  erfc  Y  .  C-l 

2  pt 

Hence,  by  substituting  C-l  into  the  Y-momentum  equation, 
integrating,  and  satisfying  the  boundary  condition  on  pres¬ 
sure  at  infinity  gives 


p  as  -b (sin  x)  0gp  (8t)  (iCt )  i3  erfc  — — —  C-2 

2  hct 


It  is  noticed  that  the  pressure  does  not  vary  with  X  (ioe,, 


A 


.  l>  .. 


there  is  no  variation  of  the  departure  of  pressure  from 
hydrostatic  conditions  in  the  lateral  direction  at  a  given 
time) . 
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The  velocity  distribution  may  now  be  obtained  from  the 
X-momentum  equation  by  firstly  substituting  C-l  and  C-2, 
resulting  in 

2 

lr  =  (COS  x)  PBg (4t)  i2  erfc  I_  +  v  .  C-3 

dt  2  /Jot  5Y2 

The  Laplace  transform  of  C-3  with  respect  to  t  after  apply¬ 
ing  the  initial  condition  U(X,Y)t_0  =  0  yields 


C-4 


sU 


the  solution  of  which  must  be  consistent  with 

U (o,s)  =  U (oo  ,s)  =  0. 

Assuming  the  particular  integral  of  C-4  is  of  the  form 


substitution  gives 


U 


P 


B<?B  (cos  y)  e 
(l-a) 


s 


3 


' 


.0  =  (a,oo)U  =  (.1,0)0 
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Hence,  seeking  the  total  solution  in  the  form 


U  =  Ae 


/f 


-M 


Y  -  Y 

+  Ee  h  e__ 

(1-a)  3 


C-5 


U(oo,s)  =  0  requires  that  A  =  0  and 

U(o,s)  ■  o  requires 


.  -BqB (cos  y) 
(1-ct)  s3 


Upon  substitution  for  the  coefficients,  the  inverse  of  C-5 
may  be  written  as 


U 


gqB (cos  Y ) (4t) 2 

(1-a) 


erfc 


erfc 


C-6 


Examining  C-6  it  is  evident  that  this  velocity  becomes  very 
small  for  values  of  x  approaching  tt/2  and  in  comparison  to 
the  velocity  of  the  main  flow,  namely 


U  = 


/c 


(1-a) 


(4t) 


5/2  r 


erfc 


2  Jict 


-  i  erfc 


2  /Jet" 


//a 


it  may  be  neglected.  Hence,  it  would  seem  reasonable  that 
the  inaccuracy  of  the  level  of  the  plate  transversely  is  not 
a  major  influencing  factor. 
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APPENDIX  D 


DATA  REDUCTION 


Considering  a  test  in  water  (as  given  in  figures  4.1, 

4.2,  and  4.3)  a  typical  procedure  used  in  reducing  the  ex¬ 
perimental  data  is  as  follows: 

Firstly,  all  strip  chart  recordings  were  con¬ 
verted  to  temperature  versus  time  curves  by  using 
figure  D.l.  Because  of  the  linear  nature  of  the 
temperature  distribution  close  to  the  surface  (see 
figure  2.2),  it  was  possible  to  determine  the  actual 
surface  temperature  at  any  time  by  extrapolation  on 
the  temperature  versus  position  (i.e.  with  respect 
to  the  surface)  curves.  Plotting  this  surface  temp¬ 
erature  distribution  on  the  temperature  versus  time 
graph  and  extending  it  back  to  the  abscissa,  permitted 
a  temporal  origin  to  be  defined.  A  typical  re-plot 
with  respect  to  this  origin  is  then  as  given  in  figure 
4.1.  Since  the  probe  was  located  at  X  =  12  +  1/8  inches, 
c  was  found  to  be  0. 34 °F/ft. sec.  from  the  slope  of 
the  surface  temperature  distribution.  Hence  knowing  c, 
the  pertinent  theoretical  curves  could  be  obtained 
but  not  before  converting 

0  =  cX  t  <t>  and  <t>  =  xt§[,  (n=l) 

o  o 
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FIG.  D.l  THERMOCOUPLE  CONVERSION  CHART 


. 


30 


by  utilizing  X  =  xXq  and  Y  =  yY  ,  to  give 


0 


cXt$ 


Similarly  it  may  be  shown  that 


1/2 


Y/2  (/ct) ly/2 


—6  o 

Selecting  /c  =  1.56  x  10  +2.0%  ft  /sec.  as  obtained 


from  tables  (19)  at  an  average  surface  temperature  and  us¬ 
ing  tables  for  $  (20)  /  a  typical  calculation  is: 


Y  =  0.045  inches,  T]  =  0.40  and  $  =  0.37 


Given 


14.1  +  2.0%  sec 


t 


0  =  0.34(1  +  1.0%)  (14.1  +  2.0%)  0.37 

=  1 . 7 9  °F  +  3.0%. 

Repeating  the  above  procedure  for  several  values  of  T)  at 
the  same  Y  and  then  for  all  Y,  results  in  the  theoretical 
curves  as  given  on  figure  4.1.  Plotting  these  curves  on 
a  temperature  versus  position  basis  resulted  in  figure  4.2. 
All  experimental  points  were  then  cross  plotted  from 


figures  4.1  and  4.2  to  give  the  non-dimensional  temperature 
profile  of  figure  4.3.  Since  overall  errors  are  best 
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revealed  on  this  type  of  curve,  all  experimental  error 
calculations  were  restricted  to  this  plot. 

Considering  a  point  in  the  maximum  velocity  region, 

0.045"  +  0.001 
12  sec  +  0.25 
1  ft.  +  0.01 
1.6  °F  +0.1  , 

0.39  +  9.0% 

0.44  +  3.0%  . 

These  percentage  errors  were  calculated  for  a  number  of 
points  in  the  region  near  the  surface,  in  the  velocity 
maximum  region,  and  also  near  the  edge  of  the  boundary 
layer.  The  maximum  errors  in  each  region  are  given  in 
the  main  text  of  Chapter  IV. 
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